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Abstract. In this summary of the application of Dyson-Schwinger equations to the theory and phenomenol- 
ogy of hadrons, some deductions following from a nonperturbative, symmetry-preserving truncation are 
highlighted, notable amongst which are results for pseudoscalar mesons. We also describe inferences from 
the gap equation relating to the radius of convergence of a chiral expansion, applications to heavy-light and 
heavy-heavy mesons, and quantitative estimates of the contribution of quark orbital angular momentum 
in pseudoscalar mesons; and recapitulate upon studies of nucleon electromagnetic form factors. 
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1. Introduction 



The world's experimental hadron physics facilities arc pro- 
viding data of unprecedented accuracy with an enormous 
potential to impact on our understanding of the basic fea- 
tures of the strong interaction; e.g., [1,2,3,4]. The expla- 
nation of that data in terms of QCD's key elements is an 
important task for contemporary theory. It is likewise vital 
for theory to link current observations to a future path of 
discovery. Success with that requires flexible tools, which 
can rapidly provide an intuitive understanding of informa- 
tion in hand and simultaneously anticipate its likely con- 
sequences. Models, parametrisations and truncations of 
QCD play this role. Such theory can also identify and ex- 
plore novel possibilities within hadron physics, whose ex- 
perimental verification could test the foundations of QCD. 

Prominent among such tools are QCD's Dyson-Schwin- 
ger equations (DSEs), truncations thereof, and models 
based on this complex of integral equations. Glancing back 
at the last proceedings volume produced for this series of 
conferences [5], it is apparent that considerable progress 
has been made in the interim. We will provide a snapshot 
of that herein, which is augmented by other contributions 
to this volume. 



* CDR is grateful to the organisers for their hospitality 
and acknowledges fruitful conversations with O. Lakhina and 
P. C. Tandy. This work was supported by: US Department 
of Energy, Office of Nuclear Physics, contract nos. DE-AC02- 
06CH11357 and DE-FG02-00ER41135; Austrian Science Fund 
FWF, Schrodinger-Ruckkehrstipendium R50-N08; and bene- 
fited from the facilities of the ANL Computing Resource Cen- 
ter and Pittsburgh's NSF Terascale Computing System. 



2. Gap equation 

In the presence of dynamical chiral symmetry breaking 
(DCSB); i.e., the generation of mass from nothing, a study 
of hadrons must begin with QCD's gap equation: 



Sip)- 1 = Z 2 (i 1 -p + 



S(p) = Z 1 J g 2 D^{p~ q )^ ltl S{ q )r:{q,p), (2) 



where J represents a Poincarc invariant regularisation 

of the integral, with A the regularisation mass-scale [6], 
D^ v {k) is the dressed-gluon propagator, r u (q,p) is the 
dresscd-quark-gluon vertex, and m hnl is the quark's A- 
dependent bare current-mass. The vertex and quark wave- 
function renormalisation constants, ifi^C 2 , -A 2 ), depend 
on the gauge parameter. 

The solution of the gap equation can be written 

S(p) = Z(p 2 X 2 )/l*l-P + M(p 2 )} (3) 

wherein the mass function, M(p 2 ), is momentum-dependent 
but independent of the renormalisation point. The solu- 
tion is obtained from (1) augmented by the renormalisa- 
tion condition 

S(P)-\, =C =ij-p + m(( 2 ), (4) 

where m{C, 2 ) is the renormalised (running) mass: 

Z 2 (( 2 , A 2 ) m bra (T) = Z 4 (C 2 , A 2 ) m(C 2 ) , (5) 

with Z^ the Lagrangian-mass renormalisation constant. In 
QCD the chiral limit is strictly defined by [6] 



Z 2 (C 2 ,A 2 )m bm (A) = 0,WA A -»C 2 



(6) 



2 



M. S. Bhagwat et al.: Mind the gap 



which states that the renormalisation-point-invariant cur- 
rent-quark mass rh = 0. 

In connection with (1) & (2) it has recently been estab- 
lished [7] that on a bounded, measurable domain of non- 
negative current-quark mass, realistic models for the ker- 
nel can simultaneously admit two inequivalent DCSB solu- 
tions and a solution that is unambiguously connected with 
the realisation of chiral symmetry in the Wigner mode. 
The Wigner solution and one of the DCSB solutions are 
destabilised by a current-quark mass and both disappear 
when that mass exceeds a critical value. This critical value 
also bounds the domain on which the surviving DCSB so- 
lution possesses a chiral expansion and can therefore be 
viewed as an upper limit on the domain within which a 
perturbative expansion in the current-quark mass around 
the chiral limit is uniformly valid for physical quantities. 
This critical mass is typically m cr ~ 60-70 MeV, which for 
a flavour-nonsinglet _ meson constituted of equal mass 
quarks corresponds to a mass m - ~ 0.45 GeV [8]. In 
arguing this case, properties of the two DCSB solutions 
of the gap equation that enable a valid definition of {qq) 
in the presence of a nonzero current-mass were employed. 
The behaviour of this condensate indicates that the essen- 
tially dynamical component of chiral symmetry breaking 
decreases with increasing current-quark mass, following 
the trend predicted by the constituent-quark cr-term ([9], 
Sec. 5.2.2). 



3. Symmetry preserving truncation 

One is forced to begin with (1) & (2) by the axial- vector 
Ward-Takahashi identity; viz., 

P M r 5Al (fc;F) = S(fc + )- 1 i 75 + z 75 S(fc_)- 1 

-2im(C 2 )r 5 (k;P), (7) 

written here for a quark and antiquark of equal current- 
mass, wherein the axial vector vertex is determined by the 
inhomogeneous Bcthe-Salpeter equation 

[r 5tl (k; P)] tu = z 2 [757^ + f A [x(q; P)UK£{ q , k- p), 

Jq 

(8) 

with X (k; P) = S(k+)r 5fl (k; P)S(k-), k± = k ± P/2 and 
the colour- and Dirac-matrix structure of the elements 
in the equation is denoted by the indices r,s,t,u. In (8), 
K(q, k; P) is the fully-amputated quark-antiquark scatter- 
ing kernel. If one knows the form of K then the nature of 
the interaction between quarks in QCD is completely un- 
derstood. The pseudoscalar vertex in (7) is prescribed by 
an analogue of (8). 

Every pseudoscalar meson appears as a pole contribu- 
tion to the axial- vector and pseudoscalar vertices [6]; viz., 

r 5 ,(k;P) P ^= lL r 7r „(fc;P)+ r^(k;P), 

(9) 



ir b (k- p) ~= ro ~" I r ff „ (fc; P) + or s (fc; P) , 

P z + m* n 

(10) 

where the _T re s are regular in the neighbourhood of this 
pole, r Vn (k;P) represents the bound state's canonically 
normalised Bethc-Salpeter amplitude: 

r,„ (fc; P) = 7 5 [iE %n (k; P) + 7 • PF^ (k; P) 
+ 1 -kk-PG 7Tn (t 1 P)+a^k^P u H n Jk;P)}, (11) 

and 

U n P„ = Z 2 tv[ 757mX^(<7;- p ): (12) 

Jq 

r A 

w.lC) =z 4 ti j 5 x. n (q;P), (is) 

Jq 

wherein the Bethc-Salpeter wave function is 

Xn n (k;P) = S(k + )r„ n (k;P)S(k-) (14) 
= 75 (fc; P) + 7 -PF« n (k;P) 
+ j-kk-Pg nn (k;P) + a liV k„P v n Vn {k;P)}. (15) 

In (9) - (15), ttq denotes the lowest- mass pseudoscalar and 
increasing n labels bound-states of increasing mass. 

Equation (7) is a statement of chiral symmetry and 
the manner by which it is broken, explicitly and dynami- 
cally. It relates the solution of the two-body problem, (8), 
in quantum field theory to the one-body problem, (1). 
QCD is violated in any approach that does not preserve 
(7). A weak-coupling expansion guarantees (7). However, 
one cannot study bound-states in perturbation theory, nor 
QCD's fundamental emergent phenomena; viz., confine- 
ment and DCSB. Fortunately, as related in [5], there is 
at least one symmetry-preserving truncation of the DSEs 
that is nonperturbative in the coupling [10, 11, 12, 13]. This 
fact has enabled the proof of exact results in QCD and 
importantly in addition their illustration using practical 
truncations to which the corrections are quantifiable. 

An exemplar: for a flavour-nonsinglet pseudoscalar me- 
son 1 [6,14] 

k< =2m(C 2 )p 7r „(C 2 )- (16) 

The Gcll-Mann-Oakes-Renner relation is a corollary of 
(16). Moreover, in deriving this expression, no assump- 
tions are made about the current-quark mass of the con- 
stituents. Hence an analogue is equally valid for a me- 
son containing one heavy-quark, a heavy-light system, and 
also heavy-heavy mesons. 

Furthermore, the validity of (16) is not restricted to 
the ground state. This entails [14] that in the chiral limit 

/° B =0,Vn>l; (17) 

For notational simplicity we've written this identity for 
mesons constituted from a quark and antiquark with the same 
current-mass. The unequal mass case is little different. 
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viz., Goldstonc modes are the only pseudoscalar mesons 
to possess a nonzero leptonic decay constant in the chiral 
limit when chiral symmetry is dynamically broken. The 
decay constants of all other pseudoscalar mesons on this 
trajectory e.g., radial excitations, vanish. This result is 
consistent with model studies; e.g., [15,16,17,18,19]. On 
the flip side, in the absence of DCSB, the leptonic decay 
constant of each such pseudoscalar meson vanishes in the 
chiral limit; namely, (17) is true V?i > 0. 

From the perspective of quantum mechanics, (17) is 
a surprising fact. The leptonic decay constant for S'-wave 
states is typically proportional to the wave function at the 
origin. Compared with the ground state, this is smaller for 
an excited state because the wave function is broader in 
configuration space and wave functions are normalised. 
However, it is a modest effect; e.g., consider the e + e~ de- 
cay of vector mesons, for which a calculation in relativis- 
tic quantum mechanics based on light-front dynamics [20] 
yields f Pl /f Pa ~ 0.6, consistent with the value inferred 
from experiment. [The expression for fy is analogous to 
(12), see [21].] As apparent in a recent lattice simulation 
[22], (17) can be realised in a framework if, and only if, 
(7) holds true and is veraciously realised. 

The two-photon decay of h mesons provides another 
important example. In the presence of DCSB the ground 
state neutral pseudoscalar meson decays predominantly 
into two photons. So long as one employs a nonpertur- 
bative truncation that preserves chiral symmetry and the 
pattern of its dynamical breaking, the ground state's two 
photon decay is described in the chiral limit by a coupling 
[23,24,25,26] 



:= % T o(-ml n = 0, 



f -<»-i7- 

z Jir 



(18) 



where T is the scalar function appearing in the matrix 
element. Equation (18) is the most widely known conse- 
quence of the Abelian anomaly. 

Now, given that /V^o = in the chiral limit, it is nat- 
ural to ask whether the Tr n ^o — > 77 transition is affected. 
Since rainbow-ladder is the leading order in a symmetry 
preserving truncation, it was used in [27] to provide a 
model-independent analysis of this process for arbitrary 
n. Therein, amongst other things, (18) is generalised and 
it is proved that 



%o(-ml n ,Q 2 ) 

Q 2 »£ Qm 4tt2 r/ ffn 

3 



m 7 Q7< 



Q 2 



,(19) 



where: 7 is an anomalous dimension; u> Vn is a mass-scale 
that gauges the momentum space width of the pseudoscalar 

meson; and Fn^ (—m 2 n ) is a structure-dependent constant, 
similar but unrelated to fn n - 2 It is now plain that Vn > 1 



lim T„o {—m n 

rh—>0 n n 



4 »A 



> 4tH 
3 



^-^ 2) (-<) 



(20) 



With the interaction described in [27], F x '(— m 
-(qq}°, and it is generally nonzero in the chiral limit. 



namely, in the chiral limit the leading-order power-law in 
the transition form factor for excited state pseudoscalar 
mesons is 0(1/Q 4 ). 



4. Heavy quarks 

The implications of (16) for heavy-light mesons are treated 
in [21,28,29]. Therefore herein we focus on heavy-heavy 
mesons, for which (16) is equally valid. For quark flavour 
Q, a constituent-quark spectrum-mass is defined via 



M Q {p 2 = C ), C 



•>pa 



-M 2 

4 Pa 



(21) 



where Mq(p 2 ) is the renormalisation-point-invariant dres- 
scd-quark mass function in (3) obtained as the solution of 
(1) & (2) with the Q-quark current-mass in (4). 3 As itlq 
is increased, Mq becomes equivalent to the so-called pole- 
mass in the effective field theory for quarkonium systems; 
i.e., non-relativistic-QCD (NRQCD). Now, with 



2M| [1 



-QQ 



/Mi 



(22) 



where e®® is a binding-energy that does not grow with 
M|, and using the renormalisation-group-invariance of 
m{( 2 )p$ Q (( 2 ), (16) predicts 



(23) 



This is an identity between the pseudoscalar and pseu- 
dovector projections of the meson's Bethe-Salpeter wave 
function at the origin in configuration space. Each ele- 
ment in (23) is gauge invariant and renormalisation point 
independent. The identity is exhibited in the heavy-quark 
limit of potential models for quarkonia; e.g., [31,32]. 

In heavy-light systems it can be shown algebraically 
[21,28,29] that (23) is realised via 

cx (m^) 1 / 2 & /«« oc l/(m^ 2 . (24) 

These results follow for heavy-light mesons because the in- 
tegrands in (12) & (13) can in this instance be accurately 
approximated via an expansion in e^ q /M% and wQ«/M%, 

where w® q is the width of the meson's Bethe-Salpeter 
wave function, which we define as the value of the rel- 
ative momentum whereat the first Chebyshev moment of 
the amplitude £^(k;P) falls to one-half of its maximum 
value. In the heavy-light meson k ~ w®^ is the typical 
momentum of the light-quark. Moreover, obtains a 
finite nonzero value in the limit Mq — » 00. It follows that 
a heavy-light meson is always of nonzero spatial extent. 

This is not true for heavy-heavy systems, as is apparent 
in Fig. 1, which depicts the evolution of the spatial size of 



In solving the Bethe-Salpeter equation, an amplitude 
peaked at zero relative four-momentum weights this value of 
Mq(p 2 ) most heavily in the dressed-quark propagator. In cal- 
culations based on [5], M c s = 2.00 GeV & M 6 S = 5.34 GeV. 
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Fig. 




calculated using the 



interaction model in [5]. Solid curve - described in (25). Dashed 
curve - linear fit to calculated result in neighbourhood of . 
Dashed vertical lines mark, from left, 1/Aff & 1/M C S . , 
approaches zero as M9 — > do. 



a heavy-heavy meson with constituent-quark spectrum- 
mass. 4 The curve in the figure is 



,<2Q 



7m 
M S Q 



In 



tm 



A 



QCD 



j M = 0.68 , T M = 8.56 , 



(25) 

with Aqcb = 0.234 GcV [5]. Plainly, with increasing con- 
stituent-mass a heavy-heavy system becomes "point-like" 
in configuration space and hence dclocalised in momentum 
space. The evolution with Mq of an observable such as 
jQQ ma y therefore be sensitive to the QQ interaction over 
a wide range of momentum scales and hence a useful probe 
of that interaction. 

In NRQCD the matrix elements for various spin states 
of a given quarkonium system are equal up to corrections 
of order Vq ~ /Mq) 2 , where k ~ W®® is the typical 
magnitude of the heavy-constituent's three- momentum in 
the meson's rest frame [33]. In this picture, 0~ + and 1~~ 
mesons, which differ because spins are anti-aligned in the 
pscudoscalar and aligned in the vector, become degenerate 
in the limit Mq — ► oo and their leptonic decay constants 

become identical; i.e., f®Q = f9®- It is noteworthy, how- 
ever, that (25) gives v 2 ~ 0.27 and v 2 ~ 0.18. Moreover, 
Vq falls only as a^(Mg). Hence, a quantitative discrep- 
ancy between fQQ & fQ® can conceivably persist until 
rather large quark masses. 

The evolution to mass degeneracy is exemplified in 
[13], which begins with a model for the dressed-quark- 
gluon vertex that appears in (2). Since the model's dia- 
grammatic content is explicitly enumerable, a symmetry- 
preserving dressed-quark Bethe-Salpeter kernel could be 
explicitly constructed. The study showed that with ris- 
ing current-quark mass the rainbow-ladder truncation pro- 



4 The evolution was calculated using the renormalisation- 
group-improved rainbow-ladder DSE truncation with the in- 
teraction model in [5] . All corrections to this truncation vanish 
in the heavy-heavy limit; e.g., see [13]. 




2.0 3.0 
£ m (19) [GeV] 

Fig. 2. Mam panel. Width/charge-squared ratio in (26). Cir- 
cles - result obtained in the approach of [5]; diamonds - data 
summarised in [35] and updated from [36]; horizontal dotted 
line indicates the conjecture of (26); and dotted vertical lines 
indicate, from left, the u, s, c, b current-quark masses fixed 
at a renormalisation scale £ = 19 GeV. Inserted panel. Circles 
- calculated evolution of fSP- Solid line - straight line fit to 



large-mq results. The c k 
by vertical dashed lines. 



b current-quark masses are indicated 



vides an increasingly accurate estimate of the mass of a 
heavy-heavy system, and the mass splitting between vec- 
tor and pseudoscalar meson masses vanishes. With the 
&-quark mass fitted to give mr(is) = 9-46 GeV, the model 
predicts m nb = 9.42 GcV. 

Little is known experimentally about heavy-heavy 
mesons. However, 1 __ mesons are readily produced exper- 
imentally and much studied. We therefore consider vector 
meson leptonic decays, for which it was noted [34,35] that 



rQQ 

Pa 



12.4 keV, 



(26) 



where (cq) is the mean electric charge of the valence- 
quark constituents in units of the electron's charge. (See 
Fig. 2.) From this empirically based conjecture one might 
conclude that on the experimentally accessible domain of 
current-quark masses f®Q ex (m^) 1 ^ 2 . This is a marked 
departure from the behaviour in heavy- light systems, (24), 
which is independent of the Qq interaction. Furthermore 
[35] , Coulomb-potential models typically give f®Q cx m®® 

whereas a linear potential produces fQ® ~ constant. Here 
is confirmation that the properties of quarkonia are a 
probe of the QQ interaction. 

In Fig. 2 we depict the result obtained for the ratio in 
(26) using the renormalisation-group-improved rainbow- 
ladder DSE truncation with the interaction model of [5]. 
The interaction is precisely that of QCD in the ultravio- 
let; namely, colour-Coulomb. However, a single parameter 
is employed to express a model for the long-range part 
of the quark-quark scattering kernel. That parameter is a 
gluon mass-scale m g = 720 MeV. It was chosen in order 
to fit / Wo & p Wo and hence the results in Fig. 2 are an un- 
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tuned prediction. Quantitatively, that for the c-quark is 
good. Indeed, for ground-state vector mesons up to and 
including J/tp the standard deviation between the calcu- 
lated width and experiment is < 15%. (NB. The calcu- 
lated masses are used in determining the widths.) Adding 
the Y(1S), that standard deviation rises to < 25 %. Quali- 
tatively, however, the mass-dependence obtained with this 
interaction does not support (26). Our preliminary result, 
illustrated via the inset in Fig. 2, is fQ® oc m®® for masses 
in the neighbourhood of the 6-quark and beyond. This is 
Coulomb-potential-like behaviour, which may be natural 
because the interaction employed is precisely that of QCD 
in the ultraviolet. We are continuing to examine the va- 
lidity of the hypothesis in (26). 



5. Quark orbital angular momentum 

It is noteworthy that quark orbital angular momentum 
is not a Poincare invariant. However, if absent in a par- 
ticular frame, it will inevitably appear in another frame 
related via a Poincare transformation. Nonzero quark or- 
bital angular momentum is thus a necessary outcome of 
a Poincare covariant description, which is why the Bethe- 
Salpeter wave function in (15) is a matrix-valued function 
with a rich structure. 

A pseudoscalar meson naturally possesses total spin 
J = and this is expressed in the fact that Xtt„ is an eigen- 
state of the Pauli-Lubanski operator with eigenvalue zero. 
Nonetheless, in the meson's rest frame one can straight- 
forwardly decompose the Pauli-Lubanski operator into the 
sum of two terms: one measuring the angular momentum 
of the quarks and another measuring their spin. In this 
way one can show that the terms in (15) characterised 
by £ and J- are purely L = in the rest frame, whereas 
the Q and TL terms are associated with L — 1. Thus a 
pseudoscalar meson Bethe-Salpeter wave function always 
contains both S- and P-wave components. 

In this connection, it is worth recalling that E Wo in 
(11) provides only a sublcading contribution to the ultra- 
violet behaviour of the electromagnetic pion form factor 
[24]. The leading power-law behaviour anticipated from 
perturbative QCD is produced by the amplitudes associ- 
ated with F Wo and G„ [25]. Moreover, it is an identity 
relating F„ n and G nn in the ultraviolet that ensures /„■„ 
in (12) is gauge invariant, and cutoff- and renormalisation- 
point-independent [25,37]. 

We exhibit the rest-frame angular momentum content 
of a ^ meson in the following way. A Bethe-Salpeter 
amplitude is canonically normalised. There are sixteen dis- 
tinguishable terms in the associated sum; viz., an £^ n £n n 
contribution plus an £ Wri J r nn contribution, etc. In the sum 
of the squares of these terms we associate sin 2 6 Vn with the 
nondiagonal contributions, in which case sin^^ gauges 
the role played by L = 1 components in the normalisation. 
In Fig. 3, for both the ground and first radially excited 
state, we plot the bound-state mass-dependence of 6„ n 
obtained in the renormalisation-group-improved rainbow- 
ladder truncation using the model interaction described in 



n * meson mass (GeV) 
1.2 1.6 2 2.4 




n * meson mass (GeV) 



Fig. 3. Mam frame: 6\r , which gauges the rest-frame admix- 
ture of L = 1 components in the ground-state h meson's 
Bethe-Salpeter wave function, plotted as a function of the me- 
son's mass. Inset: Solid curve - analogue, 7T1 , for the first 
excited state; dashed curve, for comparison. In the chiral 
limit the ground state is naturally massless whereas the calcu- 
lated mass of the first excited state is 1.04 GeV [27]. 

[5]. For both states, angular momentum is most significant 
in the neighbourhood of the chiral limit, and decreases 
with increasing current-quark mass. Notably, as measured 
by the angle 6 7Tn , at a given bound-state mass the admix- 
ture of L = 1 components in the first radial excitation is 
roughly 15% greater than that in the ground state. Mea- 
sured as a function of the current-quark mass, however, 
the situation is reversed: 61/60 « 0.5 in the chiral limit and 
this ratio increases steadily to w 0.9 at m<g(19) = 0.5 GeV. 
Our analysis continues. 

6. Baryons 

Despite material progress with the study of mesons, the 
challenge of baryons remains. A nuclcon appears as a pole 
in a six-point quark Green function. The pole's residue is 
proportional to the nucleoli's Faddcev amplitude, which 
is obtained from a Poincare covariant Faddeev equation 
that adds-up all possible quantum field theoretical ex- 
changes and interactions that can take place between three 
dressed-quarks. In formulating and solving this problem, 
current expertise is approximately at the level it was for 
mesons ten years ago; i.e., model building and phenomenol- 
ogy. However, we are a little ahead because a great deal 
has been learnt in applications to mesons. For example, we 
have acquired a veracious understanding of the structure 
of dressed-quarks and -gluons and therefore can straight- 
forwardly incorporate effects owing to and arising from 
the strong infrared modification of the momentum depen- 
dence of these propagators. (See, e.g., Sect. 5.1 of [9].) 

A tractable treatment of the Faddcev equation requires 
a truncation. One is based [38] on the observation that an 
interaction which describes colour-singlet mesons also gen- 
erates quark-quark (diquark) correlations in the colour-3 
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(antitriplet) channel [39]. The dominant correlations for 
ground state octet and decuplet baryons are + and 1 + 
diquarks. This can be understood on the grounds that: 
the associated mass-scales are smaller than the baryons' 
masses [40,41], with models giving (in GeV) rri[ ud ] o+ = 
0.74 - 0.82, m (uu)i+ = m {ud)i+ = ™ ((M ) lt = 0.95 - 1.02; 
the electromagnetic size of these correlations is less than 
that of the proton [42] - r[ ue j] + ps 0.7 fm, from which we 
estimate r( ud ) , ~ 0.8 fm based on the p-meson/7r- meson 
radius-ratio [43,44]; and the positive parity of the cor- 
relations matches that of the baryons. Both + and 1 + 
diquarks provide attraction in the Faddeev equation. 

The truncation of the Faddeev equation's kernel is 
completed by specifying that the quarks are dressed, with 
two of the three dressed-quarks correlated always 
colour-3 diquark. Binding is then effected by the iterated 
exchange of roles between the bystander and diquark- 
participant quarks. This ensures that the Faddeev am- 
plitude exhibits the correct symmetry properties under 
fcrmion interchange. A Ward-Takahashi-idcntity-preser- 
ving electromagnetic current for the baryon thus consti- 
tuted is subsequently derived [45] . It depends on the elec- 
tromagnetic properties of the axial-vector diquark corre- 
lation: its magnetic and quadrupole moments; and the 
strength of clcctromagnctically induced axial-vector <-> 
scalar diquark transitions. 

A Faddeev equation study of the nucleon's mass and 
the effect on this of a pscudoscalar meson cloud are de- 
tailed in [46]. The lessons learnt are employed in a se- 
ries of studies of nucleon properties; e.g., the nucleons' 
cr-term in [47,48], and nucleon form factors in [49,50, 
51,52]. Of particular contemporary interest is the ratio 
H P G P E (Q 2 )/G P M {Q 2 ) [1,4]. This passes through zero at Q 2 s 
6.5 GeV 2 [50]. The analogous ratio for the neutron is pre- 
sented in [51]. In the neighbourhood of Q 2 = 0, 

Ge(Q 2 ) _ _ r An 2 (27) 

where r n is the neutron's electric radius. The calculation 
shows this to be a good approximation for r^Q 2 < 1, with 
which the data [53] are consistent. It is notable that, just 
as for the proton, the small Q 2 behaviour of this ratio is 
materially affected by the neutron's pion cloud. 

Pscudoscalar mesons are not pointlikc and therefore 
their contributions to form factors diminish in magnitude 
with increasing Q 2 . It follows therefore that the evolution 
of H P G P E (Q 2 )/G P M (Q 2 ) and [i n G E (Q 2 )/G n M (Q 2 ) on Q 2 > 
2 GeV 2 are both primarily determined by the quark-core 
of the nucleon. While the proton ratio decreases uniformly 
on this domain [49,50], [51] predicts that the neutron ratio 
increases steadily until Q 2 ~ 8 GeV 2 . 

As with mesons, Sect. 5., in a Poincare covariant treat- 
ment the nucleon's quark core is necessarily described by 
a Faddeev amplitude with nonzero quark orbital angular 
momentum. The Faddeev amplitude is therefore a matrix- 
valued function that, in a baryons' rest frame, corresponds 
to a relativistic wave function with S-wave, P-wave and 
Z)-wave components [54]. In form factor studies [49,50,51] 



there is some quantitative sensitivity to the electromag- 
netic structure of the diquarks. However, the gross features 
of the form factors are primarily governed by correlations 
expressed in the nucleon's Faddeev amplitude and, in par- 
ticular, by the amount of intrinsic quark orbital angular 
momentum [55]. The nature of the kernel in the Faddeev 
equation specifics just how much quark orbital angular 
momentum is present in a baryon's rest frame. 

We see a baryon as composed primarily of a quark core, 
constituted of confined quark and confined diquark corre- 
lations, but augmented by 0~ meson cloud contributions 
that are sensed by long wavelength probes. Short wave- 
length probes pierce the cloud, and expose spin-isospin 
correlations and quark orbital angular momentum within 
the baryon. The veracity of this description makes plain 
that a picture of baryons as a bag of three constituent- 
quarks in relative S'-waves is profoundly misleading. 

7. Prospect 

Two emergent phenomena are primarily responsible for 
the observed properties of hadrons: confinement and dy- 
namical chiral symmetry breaking (DCSB). They can be 
viewed as an essential consequence of the presence and 
role of particlc-antiparticle pairs in an asymptotically free 
theory and therefore can only be veraciously understood 
in relativistic quantum field theory. The Dyson-Schwingcr 
equations (DSEs) provide a natural framework for the ex- 
ploration of these phenomena. The DSEs are a generating 
tool for perturbation theory and thus give a clean con- 
nection with processes that are well understood. More- 
over, they admit a systematic, symmetry preserving and 
nonperturbative truncation, and thereby give access to 
strong QCD in the continuum. On top of this, quantitative 
comparisons and feedback between DSE and lattice-QCD 
studies are today proving fruitful. 

Dynamical chiral symmetry breaking (DCSB) is a sin- 
gularly effective mass generating mechanism. It is under- 
stood via QCD's gap equation, the solution of which deliv- 
ers a quark mass function with a momentum-dependence 
that connects the pcrturbative and nonperturbative, con- 
stituent-quark domains. Despite the fact that light-quarks 
are made heavy, the mass of the pscudoscalar mesons re- 
mains unnaturally small. That, too, owes to DCSB, ex- 
pressed this time in a relationship between QCD's gap 
equation and those colour singlet Bcthc-Salpctcr equa- 
tions which have a pseudoscalar projection. Goldstone's 
theorem is a natural consequence of this connection. 

The existence of a sensible truncation scheme enables 
the proof of exact results using the DSEs. That the trun- 
cation scheme is also tractable provides a means by which 
the results may be illustrated, and furthermore a practical 
tool for the prediction of obscrvablcs that are accessible 
at contemporary experimental facilities. The consequent 
opportunities for rapid feedback between experiment and 
theory brings within reach an intuitive understanding of 
nonperturbative strong interaction phenomena. 

There are indications that confinement may be ex- 
pressed in the analyticity properties of the dressed prop- 
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agators. (Sect. 4.2 in [9].) To build understanding it is es- 
sential to work toward an accurate map of the confinement 
force between light-quarks and elucidate how that evolves 
from the potential between two static quarks. Among the 
rewards are a clear connection between confinement and 
DCSB, an accounting of the distribution of mass within 
hadrons, and a realistic picture of hybrids and exotics. 
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